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Abstract 

We present a systematic study of the thermodynamics of two and three-dimensional 
generahzed Lennard- Jones {LJ) systems focusing on the relationship between the 
range of the potential, the system density and its dimension. We found that the 
existence of negative specific heats depends on these three factors and not only on 
the potential range and the density of the system as stated in recent contributions. 
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1 Introduction 



In the last few years a lot of attention has been paid to Hamiltonian systems in 
which the interaction potential range is of the order of the size of the system. 
Typical examples include nuclei, metallic clusters and galaxies. In the first two 
cases the system under consideration is small comprising just a few (or less) 
hundreds of particles, while in the third one the range of the potential is long 



BO). 



One of the main consequences of dealing with long-range Hamiltonian systems 
is the appearance of negative specific heats. Negative specific heats appear in 



the 



iterature of small systems both in experimental (p; la) and theoretical 



d?!; M) studies. In particular, in a series of recent papers, we have found such 
a behavior for highly excited small Lennard Jones aggregates fl47 particles) 
both free to expand and constrained in an spherical volume (9tll(l). In this case 
the standard (12 — 6) LJ systems in 3 dimensions qualifies as a short-ranged 
potential, since the exponent of the attractive term of the potential {a = 6) is 
greater than the dimension of the system {d = 3). On the other hand, in (jlj) it 
is claimed, based on the numerical analysis of two-dimensional systems, that a 
negative specific heat region is present only if long-range forces are present. It 
is therefore necessary to make a complete study in order to clarify the interplay 
between the potential range and the thermodynamics of the system. 



The presence of such a long range interaction in the system calls for a proper 
thermodynamical description. One can, for instance, try to use Tsallis statis- 
tics llli . However, "small" systems can be also studied via micro-canonical 
statistics (8). Our studies and conclusions will be referred to the N,V,E en- 
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semble. 



This work is organized as follows. In Section 2 we present the model, the 
interaction potential and the observables that we use in the calculations. In 
Section 3 we show the results of our simulations for two and three-dimensional 
systems. Finally in Section 4 we present our conclusions. 

2 The model 

The system under study is composed by a gas of 147 particles confined in a 
spherical box, defined by the Hamiltonian H = K + V / N + Vwaiis) where K is 
the kinetic energy and the interaction potential is given hy V = T,v{rij) with 
v{r,,) = C^[^-'' - ^-1(0 < a < 12) and C„ = e(12iVa")i/(i2-)/(i2 - a). 
The Lennard- Jones like potentials as a function of r are shown in Fig. 1 for 
d = 3, N = 147, and several values of a. We used a spherical confining wall via 
an external potential Vwaii ~ (r — r^aii)^^"^ with a cut off distance Vcut = lo". 
Energies are measured in units of the potential well (e), and the distance at 
which the potential changes sign (a), respectively. The unit of time used is 
to = ya%V48e. 

The nonextensive scaling parameter, introduced in (j^), iV = 1+d ^ drr'^~^~^ 
is convenient to make the Hamiltonian formally extensive \/a/d {d is the di- 
mensionality of the system). 

n 

The model is exactly the one introduced in (jlii), and the particular case a = 6 
recovers the celebrated Lennard- Jones model. 

The set of classical equations of motion were integrated using the velocity 
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Verlet algorithm, which preserves volume in phase space ()12il: with an inte- 
gration time step between O.OOlto and O.Olto which guaranteed a conservation 
of the energy not worse than 0.01%. Initial conditions were constructed from 
the ground states of LJ systems in r — space and rescaling velocities with a 
Maxwellian distribution of velocities to the desired value of energy. All calcu- 
lations were performed once the transient behavior was over. 
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Fig. 1. Set of generalized LJ 3D potentials. 



2.1 Discussion about the N term 



To properly understand the inclusion of the N normalization term let's re- 
member the cluster expansion in real gases. By integrating in moments we 
end up with the proportionality of the partition function and the configura- 
tion integral 



1 f 



In the Van der Waals approximation 



V{r) = oo r < ro 

exp[—pV{r)] ~ 1 — i3V{r) r > tq 



(2) 
(3) 
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It can be easily shown that 



E - Eideal = Npa^rs < V{r) > (4) 

with < V{r) >= ^Jdrr^V{r). We see that is a generahzation of this term 
taking into account the adimensionahzation and the finite number of particles. 
As the configurational integral is ill-defined for long-range potentials (decaying 
slower than d) thermodynamic quantities An like the internal energy, the 
free energy, Gibbs energy, etc. associated with systems including long-range 
potentials scale like ^) 

iV y d l-a/d ^ ' 

By normalizing the potential energy with N the Hamiltonian is turned for- 
mally extensive. 



2.2 The observables 



The main observable to be extracted from our simulations is the caloric curve 
(CC), which is defined as the functional relationship between the temperature 
of the system and its energy in terms of the density i.e. T{E,p), from which 
we define the specific heat as 

^ dT/dWv 

It is then clear that a C^, < will be obtained if and only if the CC displays 
a loop. 



5 



It has recently been proposed that first order phase transitions would be uni- 

vocally signed by the amount of fluctuations in the different subsystems in 

I I 

which the system can be subdivided fll4D . In particular, the relative kinetic 
energy fluctuation Ak is defined as: 



A 



K 



2^2 



(7) 



where is the number of particles, ax the standard deviation of the kinetic 
energy per particle and T the temperature of the svstem. Since kinetic energy 



fluctuations and the specific heat are related by ()l3| ) 



< at >, 



d 



d 
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Negative values of the specific heat should be expected whenever A k ex ceeds 



the canonical value {Ak = 1.5e for d = 3 and Ak = 1-Oe for d = 2) (jlj; 121). 



3 Results 

3.1 Caloric Curves and Kinetic Energy Fluctuations in three dimensions 

We will first focalize in the study of three-dimensional systems. In Fig. 2 we 
show the caloric curve for a system composed of 147 particles with a = 9 
(upper left panel) and for a = 1, a long-range potential (upper right panel). 
For future reference, it is interesting to notice that the case a = 9 for 3 
dimensions would correspond to a = 6 in 2 dimension because for both the 
quotient ^ = 3. It can be seen that a loop in the caloric curve is present in 
both cases only when the system density is low (notice the density differences 



6 



^ 0.8 
u 

2 0.6 

g 0.4 

0.2 

2 

1.5 

< 1 
0.5 


-0.5 0.5 -0.4 -0.2 0.2 

energy (E) energy (e) 

Fig. 2. Caloric curves and relative kinetic energy fluctuations for a = 9, p = O.Olci"'^ 
(circles), p = 0.06(7"'^ (squares) and p = 0.3a~^ (diamonds) (left panels) and a = 1, 
p = O.OOlo""'^ (triangles), p = O.Olfi^^ (circles), p = 0.06o"~'^ (squares) (right pan- 
els). 

in both panels). For a short-ranged potential the relationship between the 
caloric curves and the constraining volume has already been clarified (see (J)): 
The presence of a loop in the caloric curve can be related to the formation 
of drops in configuration space. When dealing with long-ranged potentials 
the scenario slightly changes, even though the clusters can not be completely 
isolated, we still have weakly-interacting clusters, they just need more room 
to accommodate. We will come back to this point later. 

In the lower panels of Fig. 2 we show Ak for several densities (see caption 
for details). Whereas the symbols refer to calculations of via fiuctuations 
(Eq.7) the lines show the corresponding Ak estimations using the information 
provided by the caloric curves, through Eq.8. It is immediate that fiuctuations 
are enhanced well above the canonical value for the cases where the caloric 
curve displays a loop. 

In Fig. 3 we show the dependence of the caloric curve with the range of the 
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Fig. 3. Range dependence of the caloric curves for p = O.Olfj-3, d = 3. 
interaction for a constant density (p = 0.0 Icr^^), which correspond to a low 
density for all but the extremely long-range potential [a = 1), reflected by 
the absence of a loop in the caloric corresponding curve. Also notice that the 
energy which points the entrance of the system into the vapor branch is an 
increasing function of the total energy for long-range systems {a < 3) and 
then collapses to a constant value for the short-range cases studied. 

3.2 Statistical inefficiency corrections 

The importance of the loop in the caloric curve resides in the fact that it 
is a signal which identifies the presence of a first order phase transition for 
the corresponding density and lower ones. In Ref jioi the phase diagram for 
a pure LJ system was constructed by analyzing different observables. As a 
consequence of this, the phase diagram was divided into three density regions. 
In the low density regime the CC displays a loop whereas for higher densities 
the identification of the transition line was performed by a phase-space analy- 
sis. It is therefore relevant to study the relationship between the range of the 
potential and the density for which the loop disappears which we have called 
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'disjoint density pdis'' since it acts like a separatrix between two regions of the 
phase diagram. 

Up to this point we have identified the presence of negative specific heats via 
two independent signals: The presence of the loop in the CC and abnormally 
large kinetic energy fiuctuations A^- We define pdis as the density at which the 
signals (presence of the loop including errors and relative fiuctuations greater 
than the canonical value) either both disappear at the same time or they 
become inconsistent, i.e. a loop without big relative fiuctuations or viceversa. 
It is worth to notice that only in a region in p-space very close to pdis these 
two signals can become inconsistent. 

Special care should be taken because not only we are dealing with observables 
linked to mean values but to fiuctuations (second moments), as well. In order 
to disregard erroneous results we performed a statistical inefficiency study, a 



tool which was introduced by Jacucci and Rahman (jlTI ) to estimate errors 
from correlated data series. 

The block average procedure is used to calculate the error of a quantity K 
computed from a correlated data series. In order to establish the error, the 
data are grouped into h blocks of length n^. We expect the averages of the 
data in each block no to be correlated for sufficiently long nf,. The block size 
(^5) is used to obtain uncorrelated data (the statistical inefficiency s) for a 
quantity K using: 



1™ 2n7\ — (9) 
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where (j'^{K) and 0"^(< K >b) are given by: 

1 M 

AK) = J.T.iK^- < K >,otair (10) 
i=l 

1 ^ 

a2(< >,) = - 5](< - < >toi„,)' (11) 

< K >j is the average of every block and < K >totai is the average over the 
full data set of length b * rib = M. 
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Fig. 4. In this figure we show the results of the statistical inefficiency analysis for 
a = 9, p = 0.06(7~^. Panel a) shows the caloric curve for the raw data (full line, 
circles) and data taken via a random sampling (dashed line, squares). Panel b) shows 
the relative kinetic energy fluctuations for the same situations. In Panel c) we show 
more closely the differences between the random sampling and the coarse graining 
sampling (notice the errors in the random sampling case are bigger). Panel d) shows 
the statistical inefficiency as a function of the energy. 

Having determined the value of statistical inefficiency, the simulation runs are 
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divided into blocks of size rib, and sampled either via a random sampling, 
where a single value is taken at random from each block, the data set taken 
as the set of such values; or coarse-graining, where the data set is taken as the 
set of block averages. 

In Fig. 4 we show the results of the statistical inefficiency analysis for a = 9, p = 
0.06(j~^. We can see that the caloric curve does not change appreciably, we can 
not state the existence of the loop due to the numerical uncertainties involved 
(see panel c). Even though is more sensitive to numerical correlations 
between the data (it depends on both the first and second moments) this 
signal does not change qualitatively between the raw data and the random 
samphng. It is interesting to notice the behavior of the statistical inefficiency 
as a function of the system energy: There is a clear peak in the same region 
where the caloric curve displays a loop, showing that the transition from liquid- 
like to a gas-like system is in correspondence with an increase in the correlation 
time. 
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Fig. 5. Disjoint density as a function of the potential range for d = 3. 

Using results obtained with the above-mentioned methodology, we show in 
Fig. 5 pdis as a function of a. We can see that there is a clear distinction 

between long-range and short-range systems: For long-range systems {a < 
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3) pdis is an increasing function of a (a decreasing function of the potential 
range) which suggest that the loop acts as a pointer of "surfaces", since a 
long-range interaction system needs more volume to develop a structure with 
weakly interacting clusters. For short-range systems there is a collapse to a 
Pdis = 0.06a~^ value. 

Summarizing the results of this section , we have found from our calculations 
that in the case of constrained small three-dimensional systems negative spe- 
cific heats appear for all ranges of the interaction potential, i.e. for all values 
of the parameter a for densities below an appropriate threshold for each case. 

3.3 Two-dimensional systems 

We now turn our attention to systems in two dimensions. In the upper panel 
of Fig. 6 we show the range dependence of the caloric curves for a constant 
value of the density p = O.Olcr"^ for = 100 particles. 

It is immediate that there are differences between the three-dimensional and 
the two-dimensional case. In the former we found that if the density is low 
enough, then a loop in the Caloric Curve is present for every a {pdis being an 
increasing function of a). On the other hand, for the two-dimensional case we 
found a loop for a long-range potential {a = 2), and also for a short-range 
case {a = 3, that correspond to a short-range potential in d = 2). However 
the increase of a towards bigger values makes the depth of the loop almost 
vanish, which forbids us to identify negative values of Cy due to the numerical 
uncertainties involved. Moreover, If we only take into account the range of 




was defined 




compare what happens for 
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Fig. 6. Temperature, Relative fluctuations of the kinetic energy and fluctuations of 
the kinetic energy as a function of the energy for d = 2, p = O.Olo"^^. 

tifo-dimensional and t/iree- dimensional systems, we realize that we have two 
different behaviors for the same value of the potential range. Recall from Fig. 3 
that the curve corresponding to a = 9, = 3 exhibit a clear loop; on the other 
hand from Fig. 6 the case a = 6, d = 2 presents no loop. This show that the 
existence of the loop is not only related to the range of the potential, but to 
the dimension of the system as well. 

In brief, at variance with the three-dimensional case, in two dimensions the 
system does not display negative specific heats for all ranges of the interaction 
potential but only for a < 3, partially confirming the results of (Q). 

We now turn our attention to the analysis of the behavior of the quantity Ak for 
the two-dimensional case. In the middle and lower panels of Fig. 6 we show the 
relative fluctuations of the kinetic energy Ak and the standard deviation of the 
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Fig. 7. Temperature, Relative fluctuations of the kinetic energy, potential energy and 
standard deviation of the kinetic energy per particle as a function of the energy for 
d = 2. 

kinetic energy per particle as a function of the energy for d = 2, p = O.Olcr"^. 
We see once again that is a good signature of the presence of the loop. 
On the other hand, a is essentially the same for all curves shown (having 

n 

and not having a loop). In the 3D LJ case, it has been found (jl6l ) that this 
signal changes from displaying a loop to a monotonous increasing function of 
the energy at the critical density pc- This picture is compatible with the left 
panel of Fig. 7, where we show the density dependence of the caloric curves, 
Ak, and the standard deviation of the kinetic energy per particle for a short- 
range [a = 6, left panel) and a long-range [a = 2, right panel) case. It is 
worth to mention that it is only for the a = 6 system that we can make a 
correspondence with the phase diagram of a LJ fluid. Nevertheless the right 
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panel of Fig. 7 show us that the scenario does not change quahtatively for 
long-range potentials between d — 2 and d = 3. 

4 Conclusions 

In this work we have undertaken a detailed numerical analysis of the ther- 
modynamic behavior of finite confined systems in two an three dimensions at 
fixed NVE. Our goal was to study the behavior of the specific heat in terms of 
the density and dimension of the system. We identified negative specific heats 
by searching for loops in the corresponding caloric curves, and studying the 
relative kinetic energy fluctuations of the system, finding consistent results in 
both approaches. Special emphasis was taken in the statistical treatment of 
data due to correlations present in the time series. As a consequence of the 
above mentioned analysis we have found that the presence of negative specific 
heat is a quite general feature present in confined systems in two an three 
dimensions interacting with long-ranged and short-ranged interaction poten- 
tials. Using a LJ generalized potential in which the range is given when the 
value of the parameter a is fixed we have found that for three-dimensional 
systems there is always a value of the density for which the system displays 
a negative C„. On the other hand, when the system is two-dimensional, there 
is a CKiim such that above it the size of the fluctuations does not allow us to 
identify the presence of negative values of C^. 
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